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Abstract 

We have studied the collective plasma excitations of a two-dimensional elec- 
tron gas with an arbitrary lateral charge-density modulation. The dynamics 
is formulated using a previously developed hydrodynamic theory based on the 
Thomas-Fermi-Dirac-von Weizsacker approximation. In this approach, both 
the equilibrium and dynamical properties of the periodically modulated elec- 
tron gas are treated in a consistent fashion. We pay particular attention to the 
evolution of the collective excitations as the system undergoes the transition 
from the ideal two-dimensional limit to the highly-localized one-dimensional 
limit. We also calculate the power absorption in the long- wavelength limit to 
illustrate the effect of the modulation on the modes probed by far-infrared 
(FIR) transmission spectroscopy. 

PACS: 73.20.Dx, 73.20.Mf, 73.23.-b 
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I. INTRODUCTION 



Two-dimensional electronic systems, as found in single heterojunction interfaces or in 
more complex structures such as quantum wells, have been studied intensively because of 
their fundamental interest and for their potential applicationEHi. In particular, the two- 
dimensional electron gas (2DEG) with a spatially periodic modulation of its equilibrium 
density has attracted considerable attention both theoreticallyiHi! and experimentally&^il. 
One way of achieving these systems^ is to apply a voltage between a patterned gate and the 
2DEG, thereby inducing a one-dimensional periodic modulation of the density. As the gate 
voltage is increased, electrons tend to 'pool' near the minima of the modulating potential 
and, in the limit of strong modulation, the system transforms into a periodic array of ID 
quantum wires. It is clear that the restricted motion of the electrons along the wires will 
have a profound effect on the character of the plasma excitations in the system. Our main 
focus in this paper will be to study the evolution of these collective modes with increasing 
modulation of the 2DEG, and in particular, the crossover from 2D to ID behaviour. 

The theoretical approaches commonly used to discuss collective excitations in charge- 
density- modulated systems can be classified broadly as being either hydrodynamic in 
natureQHi2B0, or based on the random phase approximation (RPA)ll0Hl!ll3Hl!. The hydro- 
dynamic approaches are appealing because of their relative mathematical and computational 
simplicity. Although they vary in their level of sophistication, they all have in common the 
objective of describing the dynamics of the electronic system in terms of a closed set of 
equations for the density and velocity fields. At the simplest level of approximationi00, 
the current density in the electron fluid is determined by a local conductivity, proportional 
to the equilibrium electron density. Poisson's equation is then used to relate the electric field 
driving the current to the fluctuating electron density. Rather than explicitly determining 
the equilibrium properties of the electronic system, the equilibrium density is simply chosen 
to have some physically reasonable form. We shall refer to such theories as 'classical' in that 
no explicit quantum mechanical aspects are included in either the equilibrium or dynamical 
properties. 

Although the classical treatment is useful for obtaining a qualitative understanding of the 
collective modes in a modulated 2DEG, it has the shortcoming of not being able to account 
for nonlocal effects (e.g. in the magneto-conductivity tensor) which become important with 
increasing magnitude of the plasmon wavevector. Nonlocal effects in the plasmon disper- 
sion can be included in the hydrodynamic theory through the introduction of an electronic 
compressibilityl'0. In this way some information about the quantum mechanical equation 
of state can be built into the theory. In this respect, these theories can be referred to as 
'semiclassical'. The RPA on the other hand, provides a fully quantum mechanical descrip- 
tion of both the equilibrium and dynamical behaviour. However, being more general, the 
computational demands for its implementation are significantly greater. Furthermore, since 
it includes both single-particle and collective aspects, it is sometimes difficult to cleanly iden- 
tify excitations which are predominantly collective in natureS For both of these reasons, it 
is still worthwhile having available hydrodynamic based theories. 

In this paper, we investigate the collective response of a modulated 2DEG using the 
Thomas- Fermi-Dirac-von Weizsacker (TFDW) hydrodynamics previously developed to treat 
magnetoplasma excitations in three-dimensional parabolic wellsS'i! and electron ringsS. 
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One of the primary virtues of this approach is that it is based on a reasonably accurate 
description of the ground state properties of the electronic system, which are determined 
self-consistently from the minimization of the TFDW energy functional. Thus, unlike most 
earlier hydrodynamic treatments, the form of the equilibrium density is not chosen arbitrar- 
ily. Perturbations of the system away from the equilibrium state generate internal forces 
which drive the system back towards equilibrium. These forces are consistently included 
in the TFDW hydrodynamic equations used to describe the dynamics of the system. Al- 
though somewhat more sophisticated, the present approach nevertheless retains much of the 
mathematical simplicity of the usual hydrodynamic theories. 

Our paper is organized as follows. In Sec. [TT]we determine the equilibrium properties of a 
periodically modulated 2DEG within the TFDW approximation. Sec. [Ill] then provides the 
mathematical formalism needed to study the dynamics of the collective modes and includes 
a derivation of the power absorption which is typically measured in infrared transmission 
experiments. In Sec. [IV], we consider in detail the crossover from 2D to ID behavior in the 
periodically-modulated 2DEG, and finally in Sec. |V|, we present our concluding remarks. 



II. EQUILIBRIUM PROPERTIES 

The fabrication of mesoscopic devices often begins with a 2DEG in metal-oxide- 
semiconductor (MOS) structures or epitaxially grown heterojunctions. In these systems, 
the strong confinement of the electrons in the direction normal to the interface leads to a 
discrete spectrum of quantized subband states, while motion in the lateral dimensions is es- 
sentially free. In the limit that electronic excitations take place within only a single subband, 
the system can be idealized as strictly two-dimensional with no spatial extent in the nor- 
mal direction. In this situation, the electronic density has the form n(x, y, z) = n(x, y)S(z), 
where the planar density n(x, y) is defined by whatever additional lateral confining potential 
v ext (x, y) is imposed on the system. 

To determine the equilibrium density distribution n(x,y), we shall make use of the 
Thomas-Fermi-Dirac-von Weiszacker energy functional: 



E\n] = / dr 



Vnfr 
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Cl n 2 + C 2 l V 7V 1 -CW /2 
n{r) 

+ \jdvf rfr ^)n(r') + j drVeAvMv) (1) 

The first term in (|T|) is the Thomas-Fermi kinetic energy, the second term is the von 
Weizsacker correction to the kinetic energy, and the third term is the Dirac local exchange 
energy!!. For simplicity, we neglect any correlation contribution. The coefficients in (fl|) are 
given by (atomic units, e 2 /e = m* — h — 1, are used throughout) 

Ci = -, C 2 = ^, C7 3 = -^. (2) 
2' 8 ' 3V vr K J 

In the von Weizsacker coefficient C 2 , the parameter X w is chosen to have the value 0.25, which 
was found in other applications^ to provide the best agreement between the TFDW and full 
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density functional theory (DFT) calculations. The last two terms in Eq. (0) are the Hartree 
self-energy of the electrons and the interaction with the external potential, respectively. 

The equilibrium properties are obtained by finding the variational minimum of Eq. ([[]). 
This leads to the Euler-Lagrange equation 

8E\n] 

where the Lagrange multiplier (chemical potential) /i serves to fix the total number of elec- 
trons, N. By introducing a wave function according to the prescription n(r) = ip 2 (r), a 
straightforward calculation yields 

_^V 2 ^(r) + t; cfr (r)^(r)=^(r) , (4) 

where the effective potential is given by 

v cS (r) = 2C^ 2 (r) - ^C^(r) + 0(r) + v cxt (r) . (5) 

Here, 0(r) = / dr'n(r')/\r — r'| is the electrostatic potential arising from the electronic 
density n(r). Eq. (^) is a nonlinear equation in ^(r), and must be solved self-consistently. 
The required solution to Eq. (|) is the ground state wave function, hereafter denoted by ip , 
and the ground state energy eigenvalue defines /x. 

This formulation of the equilibrium properties of the electronic system is quite general 
in the sense that the physical situation is completely specified by the form of the external 
potential. For the case of ID modulation, the potential v ext (x) will be assumed to be a 
smooth, periodic function of x with period a. This potential can of course be represented 
as a Fourier series, but we shall consider only the simplest form 

Vatix) = -V COS ( — — ) (6) 



a 



having a single Fourier component. This choice is sufficient to investigate the important 
effects of modulation on the properties of a 2DEG, and has been applied successfully by 
other workers0@il to model laterally microstructured field-effect devices used in FIR ex- 
periments. 

The translational invariance of the system along the y-direction, together with the pe- 
riodicity of the external confining potential, implies that the solution to Eq. (|j) will only 
depend on x and will have the property that i/jo(x + a) = if)o(x). Thus, we can restrict our 
calculation to the unit cell x G [—a/2, a/2]. For potentials such as (^j) with inversion sym- 
metry about the center of the cell, the desired solutions are those satisfying the boundary 
conditions 

(±|) = • (7) 
To complete the determination of ipo(x), we must finally impose the normalization condition 

dx \i/j (x)\ 2 = n 2 D , (8) 

a/2 
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where fi2D is the average electron density in the modulated 2DEG. 

Before discussing the nature of the ground state solution to Eq. (||), we first note some 
technical points involved in the calculation. As is typical of DFT schemes, the solutions 
are obtained using an iterative procedure, and issues of numerical stability and convergence 
therefore arise. The usual source of numerical instabilities is the long-ranged nature of the 
Coulomb interaction. In the course of iterating, a small deviation of the charge density 
from its true value leads to a change in the effective potential that tends to induce an 
overcompensating screening charge on the next iteration. In our particular case, this results 
in an oscillation of the charge between the central and outer regions of the unit cell, which 
inhibits the desired convergence of the iterative procedure. This problem can be remedied 
by a combination of two strategies. First, we introduce an algebraic mixing of the effective 
potential 

= a v eS [n®} + (l-a)«g, (9) 

where v^g is the effective potential on the i-th iteration, is the density generated by 
this potential and f e fr[n.w] is the effective potential that this density would lead to. The 
mixing parameter a lies between zero and one and reducing its value reduces the change 
in the potential from one iteration to the next. This by itself will eliminate the instability 
in some cases, as for example when the modulation periods are relatively small. However, 
to access the physically relevant modulation periods (100 nm < a < 1500 nm in GaAs), an 
extremely small value of a would have to be used, which results in an unacceptably slow 
rate of convergence. This problem can be overcome by means of a second strategy in which 
the long-range nature of the Coulomb potential is effectively screened. The screening is 
conveniently formulated in Fourier space using the algorithm 

0" + "(q) = ^) + (l-f)^(q), (10) 

where 0(q) and n(q) are the 2D Fourier transforms of the electrostatic potential and 2D 
electron density, respectively, Q = \J q 1 + k 2 , and k is an artificial screening parameter. One 
can check that the first term on the right hand side of (1C), in which the density acts as 
a source, leads to an exponentially decaying potential in real space. Nevertheless, at self- 
consistency, (fi( l+1 \q) = 0^(q) = </>(q) = 2nn(q)/q, which is just the Fourier transform of 
the actual Coulomb potential 4>(r). By choosing k of order unity, we were able to obtain 
convergence in all cases of interest with a > 0.005. Although these values of a are still 
rather small, they provided a rate of convergence which was significantly better than that 
achieved with the algebraic mixing scheme alone. 

The solution to Eq. (f|) is shown in Fig. [I] for a range of potential modulation amplitudes. 
It is useful to represent these amplitudes in the dimensionless form V /Ef, where Ef = |fc 2 
is the Fermi energy of the uniform 2DEG and the Fermi wavevector kf is defined by fc 2 = 
27iri2D- With an application to GaAs in mind, typical densities are in the range 10 11 - 
10 12 cm~ 2 . Using the effective mass m* = 0.067m e and static dielectric constant e = 13.0 
of GaAs, the atomic unit of length is the effective Bohr radius = 103 A, and the unit 
of energy is the effective Rydberg Ry* = e 2 /2ea* = 5.4 meV. In these units, a density of 
10 n cm~ 2 is approximately 0.1(cIq)~ 2 ; this is the value of the average density used to generate 
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the density distributions in Fig. |l[ The equilibrium densities in this figure span the range 
from weak to moderate to strong modulation, and illustrate the way in which the cell- 
boundary density decreases with increasing values of Vq. In the extreme localized limit, the 
system consists of isolated quantum wires, each of which is confined within an approximately 
harmonic potential well of curvature 2n 2 Vo/a 2 . It can be seen that the density in this limit 
is similar to the semi-circular profile found in the classical electrostatic approximation!!!, 
and only deviates significantly from this form at the edges of the quantum wire as a result 
of the quantum mechanical barrier penetration included within the TFDW approximation. 

In Fig. [2], we have constructed a 'phase' diagram in the (akf-Vo/Ef) plane which marks 
the boundary between localized and delocalized solutions of Eq. (|]). To generate these 
curves, we have used the criterion no(±§) = 0.1n 2 D- This condition is rather arbitrary, and 
a value somewhat greater than or less than 1/10 could have been used to generate a similar 
set of curves. The solid and dashed curves correspond to average 2D densities of Hid = 0.1 
and n2D = 1-0, respectively and span the range of densities of physical interest. The open 
circles represent the 'phase' boundary in the case of noninteracting electrons obtained by 
ignoring the electrostatic and exchange potentials in (|^). In this limit, one can show from 
Eq. (H) that the scaled wavefunction kfipo(x) as a function of the scaled length x/a depends 
on the density only through the two parameters akf and V /Ef. As a result, the phase 
boundary in the (akf-Vo/ Ef) plane is a universal curve which is valid for any density. It 
can be seen that this scaling property of the wavefunction does not apply once interactions 
are included. We have also indicated by the labeled points the values of the parameters 
corresponding to three of the density profiles in Fig. |I[ These were obtained for a density of 
^2D — 0.1, so that the position of the points relative to the solid curve in Fig. |2| is relevant; 
point "(b)" is associated with a strongly modulated but extended density while point "(c)" 
is well into the localized region. 

The parameter akf can be used to distinguish two regimes of interest: the rapidly varying 
regime akf <C 1 where the modulation period is small compared to the Fermi wavelength 
of the average density, and the slowly varying regime akf ^> 1. In the former limit, all 
of the curves approach an asymptote which is weakly dependent on density. This limit 
of strong quantum confinement can be understood simply in terms of the criterion that 
the potential barrier height, Vq, is large compared to the quantum mechanical zero point 
energy h 2 /ma 2 . We thus conclude that the 'phase' boundary is given approximately by 
Vo/Ef ~ (akf)~ 2 , which is consistent with the numerical calculations. In this limit, the 
quantum kinetic energy is deciding the question of localization. In the opposite limit of a 
slowly varying potential (akf ^> 1), the von Weizsacker kinetic energy becomes negligible 
and we recover the Thomas-Fermi (TF) approximation. The experimental situations of 
interest are typically in this TF limit. The 'phase' boundary now approaches a straight line 
with zero slope in the noninteracting case (open circles) and finite slope in the interacting 
case (solid and dashed lines). The behavior of the noninteracting curve follows from the TF 

density, utf(x) = fi2D + cos (^p0 ; an d our criterion for localization implies Vo/Ef ~ 1, 
which is to be expected when the energy of the gas is exclusively kinetic. On the other hand, 
the behavior of the interacting 'phase' boundaries can be understood in electrostatic terms. 
The external potential v ext (x) can be viewed as arising from a modulation of the external 
positive background density about the average value n 2 D- the amplitude of this density 
modulation, An ext , is related to Vq by Vq = ^-An cxt = aAn cxt . Since the electrons simply 
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neutralize the positive background locally in the TF limit, the criterion for localization 
becomes V ~ an 2 D, which implies V /Ef ~ —\—(akf). Thus, the 'phase' boundary has a 

slope inversely proportional to the square root of the average density, consistent with the 
results in Fig. 0. 

We may also consider general solutions of Eq. (|4]) which for the case of one-dimensional 
periodic modulation will take the form of Bloch-like states in the x direction, Vw( x ) = 
e lQxX u nqx (x), where u nqx (x) is a periodic function of x, multiplied by a plane wave factor, e tqyy , 
in the y direction. Here, q x is restricted to the first Brioullin zone, — n/a < q x < ir/a, and 
n is a band index. The q = state in the lowest band is the ground state ipo(x) considered 
previously. Although the significance of these general solutions is not immediately apparent, 
we shall see that they have some relevance to the calculation of the collective excitations 
in a modulated 2DEG, which similarly exhibit a Bloch-like structure. For the purpose of 
comparison with the plasmon bands to be calculated, we show in Fig. |3] the TFDW energy 
bands in the x direction for the case of moderate modulation ("(b)" in Fig. 1). 



III. COLLECTIVE EXCITATIONS 

A. Hydrodynamic Equations 

In order to determine the plasma modes in the modulated 2DEG, we adopt the TFDW 
hydrodynamic approach developed previously^. This is based on the usual continuity equa- 
tion 

^ + V>v) = 0, (11) 

and the momentum equation 

= nF — nw x uj c , (12) 
where F = F int + F ext includes both the internal force acting on the electrons, 

F int (r,t) = -V 

as well as any additional time-dependent external force, F ext (r,t). The potential term tyj 
in Eq. (|13"D contributes the expected force corresponding to the internal TF pressure and 
Coulomb-derived potentials, while the remaining term is associated with the von Weizsacker 
kinetic energy. For completeness, we have also included in Eq. (0) the magnetic force due to 
an externally applied magnetic field, B, which we shall take to be uniform and perpendicular 
to the 2DEG (i.e., B = Bz). This force is expressed in terms of the cyclotron frequency 
vector uj c = eR/m*c. Although we defer an explicit treatment of magnetoplasma modes to 
a future paper, we shall develop the equations for this more general situation as very little 
additional effort is needed to do so. 

Self-consistent solutions to our hydrodynamic equations are obtained by linearizing the 
density in small deviations from equilibrium, viz., n — > n + 5n. If the time-dependent 



n 



dv 

Iff 



— + v ■ Vv 



v e ff(r,t) 



2 VM) 



(13) 
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density is represented as n = ip 2 , the fluctuation of the von Weizsacker wave function, Sift, is 
related to the density fluctuation by 5n = 2ip Sip. Retaining only first order quantities (the 
velocity v is itself first order), Eqs. (|11|) and fll2"|) yield 



d5n 



+ V • (nov) 







and 



— = 6F — v x uj ( 

dt 



where the fluctuating force is given by 
5F = -V 

with 

5v eS 



dv Pff - —V dip + — dip 



2ip 



2 M 



+ 5F- 



ext. 



(14) 



(15) 



(16) 



(17) 



The bracketed term in Eq. (^) will be denoted by /, so that 5F mt 

For the purpose of determining the normal mode frequencies of the system and the 
associated mode densities, 5F ext can be set to zero. We shall later consider the response 
of the system to external fields in the calculation of the power absorption. Due to the 
translational invariance in the ^-direction, all of the fluctuating variables (Sn, 5ip, f and v) 
will have the form of a propagating wave e 1 ^^ - ^' with a;-dependent amplitudes. Making 
use of this dependence, Eqs. ( O) and ( 15|) can be expressed in the form 



- iudn + in q y v y + ^-(n Q v x ) 



and 



zujSF - {uj c x SF) , 



(18) 



(19) 



respectively. In these equations, only the x-dependent amplitudes are displayed. Recalling 
that 5F = —V/, we have 



uj c x 5F = (iq y u} c f)Z 



df 

dx 



(20) 



The use of Eq. 



in Eq. (|T9"D , along with Eq. flUf) , then leads to 

dn 



— uJ c )5n = q un f — uj— \ n 



d_ 

dx 



df 
dx 



(21) 



Eq. can also be expressed in terms of the fluctuating wave function by using uq = ipQ 
and Sn = 2tpo5ip. With this substitution, we find 
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u(u 2 - u 2 c )5ip = -qluipof - uJ^'of - -en/>of - u c q y ij)' Q f , (22) 

where derivatives with respect to x are now denoted by primes. Noting that ipof" = {ipof)"~ 
f'f - 2%f, Eq. @ yields, 

2 2 ipo 

- L(VWT - ^AvWO . (23) 

The advantage of this form is that the function (ipof) and its derivatives now appear on the 
right hand side, where 

M = ipoSv eS + y (f + q y) ^ - T 6r 

= M5il) + hdifj . (24) 
The operator M is defined by MSip = ^o^eff and 

is just the Hamiltonian (for q y = 0) determining the ground state von Weizsacker wave 
function ipo(x). Substituting Eq. (|2~4]) into ( j2"3"|) yields 

u{u 2 - LU 2 )5ip = uX^hCh + M)6i/) - uj c q y ^{h + M)5tp . (26) 

In the limit of vanishing magnetic field, this reduces to 

X w u 2 5%l) = h(h + M)5ip . (27) 

Due to the periodicity induced by the modulating potential along the x-direction, the 
fluctuating part of the wave function will have the Bloch-like form 

5^ = e ^ x Y, c gVg , (28) 

G 

where c G is a Fourier expansion coefficient and tp G = ~/a e%Gx • These basis functions have 
been chosen to satisfy the orthonormality condition 



ra/2 

/ dxip* G (x)ip G ,(x) = 8 GG > , (29) 

J-a/2 



where G = ^ {n = 0,±1,±2,...) is a one-dimensional reciprocal lattice vector in the 
x-direction. Substituting the Fourier expansion ( [28] ) into fl2"6]), we obtain the equation 

u{u 2 - u 2 )c G = -uj c q y A GGI M G r G ,,C G n + U ^ B GG > M G t G n C G " , (30) 

G'G" G'G" 
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where the various matrices appearing in this equation are defined as 



A GG , = —(G - G')\mp [G - G') , (31) 



Bgg> = ^ w h G G 

1 



and 



(q x + G) 2 + q 2 6 GG > + ^^(^eff -H)[G- G'] , (32) 



M GG > = M GG > + KB GG > ■ (33) 



In obtaining ([H]) we have used the identity ^ = -^(lnifj Q ), and denote the Fourier transform 

of a periodic function with an overline, e.g. ipo(x) = J^ci'olG^G- 

The remaining matrix Mqg' i n (|33|) consists of several terms. From Eq. (p4|), we have 

MSip = ip 8v cS 

= 4C^ 2 6ij - ^C 3 ^ + ^8(j> . (34) 

The last term in ( |34"D involves the x-dependent part of the total electric potential fluctuation, 
5<p. In terms of the density fluctuation 5n, we have 



5cj)(r) = J dt 
dr 



, 8n{r') 



r-f 

, 2 Mx'y (qxX ' +qyy,) J2g> cg>Vg>{x') 
| r — r' | 

Att p iGx 

= e ^ +gy y) E 6 MG ~ G'\c GI , (35) 

a GG' yJ{q x + G) 2 + q 2 

where we have noted that the 2D Fourier transform of r~ l is (2iv/q). It is clear from this 
expression that 5<fi(r) is proportional to e lc,yV and that the x-dependent part has the expected 
Bloch-like form. With this result, (|34]) yields 



Mgc = M« G , + Mq G , + Mq G , , (36) 

where 

M gg> = -4^o[G - G'\ , (37) 



M, 



GG' 



3C,— 



2^ 



MG-G'} , 



(38) 



M. 



H 

GG 1 



An 4j [G - G"] l ip [G" - G') 



7^ 



(q x + G") 2 + q 2 



(39) 
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The superscripts K, X and H refer respectively to the kinetic, exchange and Hartree terms 
in il) Q 5veg. It is clear that the kinetic and exchange matrices depend on G — G', but the 
Hartree matrix, which is associated with a nonlocal operator in position space, has the 
property M§ +GoG , +Go (q x ) = Mg G ,(q x + G ). 

Eq. ( POD is a nonlinear eigenvalue problem which must be solved numerically for the 
eigenvalue u and eigenvector c. In principle, the dimension of the eigenvalue problem is 
infinite, and for practical purposes, a truncation of the expansion to some finite number 
of G vectors is required. However, one can always check that the results for the modes of 
interest have converged by systematically increasing the number of G vectors. Once the 
eigenvalues ui have been determined, the corresponding eigenvectors c can be substituted 
into Eq. ( pq) to determine the mode densities. 

To close this section, we note some general consequences of the equations we have ob- 
tained. In the q y — > limit, (|30|) reduces to 

(U 2 - L01)C G = B GG'M GIG nC G n , (40) 

G'G" 

and the magnetic field appears explicitly only on the left hand side of the equation. This 
implies that the magnetoplasma frequencies have the property 

f n(&, % = °! B ) = u n(<lx, q y = 0;B = 0)+ u 2 c , (41) 

that is, a simple cyclotron shift of the zero-field frequencies, u n (q x ,q y = 0; B = 0). In 
addition, the mode densities in this limit have exactly the same spatial distribution as the 
corresponding B = mode densities. 

Another general property can be deduced from Eq. fl27D by making use of the Bloch- 
state basis, ip nqx (x), introduced at the end of Sec. II. These states are eigenstates of h 

with eigenvalues e n (q) = E n (q x ) + ^ w Qy ~ A*, where E n (q x ) are the ID band energies and 
the chemical potential /i is equal to E (0), the q = energy of the lowest band (n = 0). 
Expanding Sip as 

Sip = d n ip n g x , (42) 

n 

and substituting this expansion into Eq. fl27D , we obtain the eigenvalue problem in the 
alternate form 

(q)J2l £ n{c05 nn , + M nn/ ]d n ,. (43) 

n' 

By its definition, £o(0) = for the lowest band, and as a result, Eq. ( f43"D will have a 
nontrivial solution at q = with frequency u = 0. In other words, the lowest plasmon band 
disperses from zero at q = and, according to Eq. (f4lf) , this implies that there will be a 
magnetoplasmon branch which disperses from u c . 

B. Power Absorption 

The main method for studying the collective modes in these systems is by means of FIR 
absorption experimentsE^TEaEa. To make contact with these experiments, we consider in this 
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section the calculation of the power absorption which, in addition to the mode frequencies, 
contains information about the oscillator strengths of the observed excitations. 
The instantaneous power absorption is given by 



P(t) = jrfrf d (r,t).E cxt (r,i), (44) 

where j ind (r,t) is the current induced by an external electromagnetic field. For a uniform 
radiation field polarized in the x-direction and incident normally on the sample, the external 
field is spatially uniform and is given by E ext (r,t) = ^E (e~ luJt + e lw *)x. The physically 
relevant quantity is the time-averaged power absorption which is given by 

{P)t = \E, J drRef» d (r,oj) . (45) 

Since the induced current j x nd for the situation of interest is a periodic function of x and 
independent of y, the time-averaged power absorption per unit area is given by 

(P)t 



A 

where 

1 W2 



E o Bej% A [G = 0,u] , (46) 



£*[G, = — dx e-^j^(x, uj) (47) 

Q, J— a/2 



is the Fourier coefficient of the induced current. 

The current density is determined by j x nd (x,uj) = — n (x)v x (x, uj) where the velocity is 
the solution of 

— = -TV + 5F - v x uj c . (48) 

We have included in this equation a phenomenological relaxation rate, which accounts for 
momentum nonconserving scattering processes of the electrons. Because of the relaxation 
rate, the frequency uj appearing in Eq. (fl9|) is now replaced by u = uj + ij. Noting that 5F 
now includes the additional term 5F£ xt = —Eq, we obtain the following expression for the 
current density at q = 0, 

f: d (x,uj) = -n (x)v x (x,uj) = ~ a^ rco f |^ + ff>) • ( 49 ) 

Taking the q = limit of the continuity equation in (0), along with Eq. ([E3), we find 

u{ U - u c )5n = -u- [n - j - ujE q — . (50) 

Thus, the net effect of including the external driving field E ext is to convert the eigenvalue 
problem in Eq. (|30D , into the set of inhomogeneous equations 

uj(uj 2 - uj 2 c )c g -uj B GG ,M G , G „c G » = -iujG^) [G]EQ . (51) 

G'G" 
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Substituting Eq. into Eq. (|47|), we obtain 



—2d) x „ iu;i?o 



ji nd [G = 0, u] = £ G MG)fa + ^r^MG = 0] , (52) 

v4^ 2 - uj 2 ) ^ uj 2 - uj 2 

where, in deriving this result, we have used the fact that (n f) is periodic, n = ipQ, and 
ipo[—G] = i/)q[G]. The quantity fa is the Fourier coefficient of (ipof) defined in Eq. Q23D, and 
is related to cq by the equation 

fG = Y.MGG>CG> ■ (53) 
G> 

The substitution of (|52"D into (|46[) provides our final expression for the power absorption. 
Since no[(j = 0] = \faifi2D-, the last term in Eq. (|52"D contributes an absorption peak at 
the cyclotron frequency, uj c , as found for a uniform 2DEG. This peak shifts to uj = for 
B = 0, and represents the expected Drude-like absorption associated with the resistive losses 
in the 2DEG. The other term on the right hand side of Eq. ( p2|) accounts for the density 
inhomogeneity and has the effect of reducing the amount of Drude absorption. 

The cancellation of the Drude peak can be seen more clearly by writing the power 
absorption in an alternate form. Starting with the continuity equation fll~8|) with q = 



J^l nd = 0, (54) 



jind[£ = 0) a,] = — / x Sn(x) dx + v^j ind (a/2) . (55) 

\/a J-a/2 



multiplying by x and integrating over a unit cell, we find 

ra/2 
'a J-a/2 

The substitution of this expression into Eq. fl46|) then yields 

(P) 1 1 ra/2 I 

LJ1 = — u e Im - / xSn{x, uj) dx + -E Ref* d {x = a/2, uj) . (56) 
A 2 a J-a/2 2 

In the localized limit, the boundary current j^ d (x = a/2, cj) vanishes and the power absorp- 
tion is then determined by the induced dipole moment of the charge density fluctuation. In 
this case, there is no absorption at the cyclotron frequency, oj c , despite the appearance of 
the resonant denominators in Eq. ([521) . in the opposite limit of a weakly modulated system, 
the induced dipole moment will be small and the power absorption will be dominated by 
the second term in (|55|). We then recover the Drude absorption discussed previously. 



IV. THE 2D TO ID CROSSOVER 

In this section, we study the plasma modes of a 2DEG subjected to the modulating 
potential given in Eq. (|6|). Our main interest is in the evolution of these modes as a function 
of the strength of the modulation. The crossover from 2D to ID behaviour will occur in the 
vicinity of the 'phase' boundary illustrated in Fig. 1. Once the modulation is sufficiently 
strong, the electron layer separates into an array of quantum wires and we can expect very 
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different behaviour from the original 2D situation. Although the theory in Sec. Ill was 
developed with the inclusion of a magnetic field, we shall restrict our investigation to the 
zero field limit, in which case Eq. fl3"0|) reduces to 

u 2 c G = b gg'M g > G iic G ii . (57) 

G',G" 

All the mode frequencies presented in this section are based on an analysis of this equation. 



A. The Uniform 2DEG 

It is useful for the purpose of orientation to begin with the homogeneous 2DEG (t> ex t = 0) 
as treated in the TFDW approximation. The interaction matrices in Eqs. (|37|) -(|39D, as well 
as the 5-matrix in (0), are diagonal in this case and we readily obtain from ( [57] ) the plasmon 
dispersion relation 

uj 2 (q) = 2nn 2D q + y2Cin 2D - -C-sV^d) <? 2 + "jV 

= ^+(^/-^/)^ + T^ 4 - ^ 

At long wavelengths, this gives the expected 2D plasmon frequency u>2D = v / ^™2D<? which 
has the characteristic y/q dependence. At shorter wavelengths, the effects of the TF (Ci), 
exchange (C3) and von Weizsacker (X w ) energies become of increasing importance. It is inter- 
esting to note that the exchange interaction gives a negative q 2 coefficient which counteracts 
the positive dispersion coming from the TF kinetic energy. As a result, the q 2 coefficient 
goes to zero at a density having kf = 2/ V. 



B. Weak Modulation 



We now consider a weak modulating potential corresponding to the phase-point "(a)" in 
Fig. ||, with coordinates (akf = 31.71, V /Ef = 6.4). In Fig. f| we illustrate the q x - dispersion 
for q y = 0, as determined from Eq. fl5"T|). Also shown for comparison as the solid line is 
the uniform 2DEG dispersion relation given by Eq. fl5"8"|). The similarity of the dispersion 
curves for the two cases is notable. The only obvious difference is the appearance of a gap 
in the plasmon dispersion at the zone boundary, q x = ir/a, which is associated with the 
development of plasmon bands. This problem was treated theoretically by Krasheninnikov 
and Chapliki using degenerate perturbation theory. They found that the size of the gaps 
induced by the density modulation is given by A ~ (\N n \/n2D)u n , where u n is the unper- 
turbed plasma frequency for wavevector q x = q n = mc/a and N n is the nth coefficient in 
the Fourier expansion of the 2D electron density (N n = no[2irn/ 'a]/ ' ^fa in our notation; odd 
values of n give rise to zone boundary (ZB) gaps, while even values give rise to zone center 
(ZC) gaps). Since their analysis is based on a different formulation of the problem and 
includes only the Hartree interaction, we present for completeness a perturbation theory 
calculation of the gaps within the TFDW hydrodynamics in Appendix [A|. We find that the 
gaps are determined by the equation 
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c4 = ^o(9n) ± fw(j(g n ) 



3 /= 2,3. 4 

- -CV^D^n + -A„,g n 



n [27m/ a] 



(59) 



which is consistent with the result of Krasheninnikov and ChaplikB in the long wavelength 
limit where the Hartree interaction is dominant. In Fig. [15] of Appendix |A|, we show some of 
the Fourier coefficients n^[27m/a] as a function of Vo/Ef. The n^[27r/a] coefficient is linear 
in Vq whereas the higher Fourier components are at least of order Vq and are therefore much 
smaller in the weakly modulated regime. This accounts for the relatively small magnitude 
of the gaps in Fig. f| as compared to the lowest zone boundary gap. 

In Fig. H we show the plasmon dispersion in the y direction (i.e., perpendicular to the 
modulation direction) for modes at the zone center (q x = 0, solid circles) and at the zone 
boundary (q x = ir/a, open circles). Specific modes of interest that will be referred to in 
the text are indicated by encircled numbers. Of the q x = modes, there is one which 
starts at zero frequency (mode "1"), as explained at the end of Sec. |III A| . To determine 
the small-g^ behaviour analytically, we make use of the long-wavelength form of Mqc 
derived in Appendix B. It is shown there that the 2D plasmon dispersion for q x = with 
q y — > is given by u 2 ~ 27in 2 Dq y , the usual plasma frequency of a uniform 2DEG. This 
is a general result independent of the degree of modulation of the 2DEG. In other words, 
the modulation of the 2DEG has no effect on the long-wavelength dispersion of plasmons 
in this direction. However, as discussed in more detail in Sec. |LV Q this is not the case for 
plasmons propagating in the x direction. 

There are several other notable features in Fig. |5|. The first concerns the modes at 
Qx = it /a which behave as 5n(x + a) = —5n(x) and have a periodicity of 2a. The two lowest 
modes (the "2-8" and "3-7" branches) indicated by open circles in Fig. | cross close to q y — 1. 
This implies that the gap at the zone boundary first decreases to zero as a function of q y , 
and then increases again. Such a crossing is possible since the two modes in question have 
different symmetries. To illustrate this more clearly, we show the induced charge density 
fluctuations for various modes in Fig. ^|. The left panel of this figure corresponds to q y = 
while the right panel is for q y = 3ir/a. The sequence from top to bottom illustrates the five 
lowest modes, starting with the q x = mode in the top panel, followed by the two q x = n/a 
modes in the middle panel and finally the next pair of q x = modes in the bottom panel. 
The right panel of Fig. |] shows how these mode densities change when q y is increased to 
37r/a. It is the pair of q x = rr/a modes in the middle panel which exhibits the frequency 
crossing in Fig. [|. For q y values up to the crossing point, the lower frequency mode (the "2-8" 
branch) is symmetric with respect to reflections about x = 0, while the next highest mode 
(the "3-7" branch) is antisymmetric. The ordering in frequency of these two modes is then 
reversed at the crossing point, with the odd-parity mode lying lower. Similar behaviour was 
found for the model densities with weak modulation considered in Ref. [ El, but for strong 
modulation an anticrossing behaviour was observed. As we shall see, the crossing behaviour 
we find persists for even stronger modulations (see Fig. |9[). 

A further examination of Fig. |5| shows the odd-parity mode at q x = ir/a approaching 
the lowest-lying q x = mode with increasing q y (points "6" and "7"). The reason for this 
can be explained by comparing the mode density ( "6" ) in the top-left panel of Fig. with 
the odd-parity density ( "7" ) in the middle-right panel. Both mode densities are seen to be 
localized at the cell boundaries where the equilibrium density is lowest, and both have a 
very similar spatial profile in this region. The near degeneracy of these modes for large q y 
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indicates that there is a weak interaction between the density fluctuations localized at adja- 
cent cell boundaries and that as a result, these modes effectively propagate independently 
of each other. This behaviour is analogous to the situation in a metallic slab where the 
symmetric and antisymmetric surface plasmons become degenerate at large wave vectors. 
The localization of the mode density at the cell boundaries increases with increasing q y , 
implying that the modes are channeled in their propagation along the low density part of 
the gas. 

We also see a similar convergence of the other q x = ir/a mode with the next higher 
q x = mode at large q y . The explanation for this is the same as given above. The density 
fluctuation "8" shown by the solid line in the right-middle panel of Fig. |6] is the even-parity 
version of the odd-parity density fluctuation "9" shown by the dashed curve in the lower- 
right panel. A similar pairing of q x = and q x = ir/a modes would also be expected for the 
higher lying modes in Fig. [|. This behaviour is in fact more evident when the modulation 
amplitude of the equilibrium density increases (see Figs. |^ and |T^). We should emphasize 
that the convergence of pairs of q x = and q x = ir/a mode frequencies implies that these 
modes exhibit a very weak dispersion with respect to q x . The absence of an effect of the 
e lQxX phase modulation from one unit cell to the next reflects the lack of interaction between 
adjacent density fluctuations. 

Returning to Fig. |], we see that the gaps at the zone center are unobservably small for 
this case of weak modulation. However, each of the higher lying modes at q x = corresponds 
to two distinct modes. This becomes apparent for the first excited q x = modes in Fig. |5| 
(starting at the point labeled "4" and "5") where it is seen that the frequencies separate with 
increasing q y . The lower-left panel in Fig. [] shows the mode densities for this pair of modes 
at q y = 0, and the lower-right panel shows the densities at q y = 3n/a. One mode has even 
parity with respect to the center of the cell, and the other has odd parity. The latter mode 
eventually evolves with increasing modulation into the lowest odd-parity mode of an isolated 
wire. Since this mode has a finite dipole moment, it will couple to an external radiation 
field and will contribute to the power absorption as indicated by Eq. (|56|). To illustrate 
this, we show in Fig. ^ the calculated power absorption for different modulations of the 
2DEG. The lowest curve labeled "(a)" corresponds to the case of weak modulation being 
considered here and shows a small peak at the frequency of the first excited q x = mode 
"4". This is essentially a bulk 2D plasmon at a wavevector q = 2ir/a. In principle, other 
q x = odd-parity modes should be observable, but their oscillator strengths are too small 
to show up in the power absorption. In fact, the dominant feature in the power absorption 
of curve "(a)" is the strong Drude peak at uj = which is to be expected since the system 
is only weakly perturbed from a homogeneous 2DEG. We shall return to a more systematic 
discussion of the power absorption later. 

C. Moderate Modulation 

As the modulation is increased, the equilibrium density profile becomes more localized 
(see the curve labeled by "(b)" in Fig. [I]). The effect of this increased modulation on the 
wavevector dispersion along the x-direction is illustrated in Fig. |^. A comparison of this 
figure with Fig. ^ reveals several notable changes, the most dramatic being the increase in 
the magnitude of the gaps at both the zone center and zone boundary. As we have already 
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discussed in Sec. [IV B| , the size of the gap is related to the magnitude of the density Fourier 
coefficients, so this result is to be expected. We also notice that the q x - dispersion of the 
lowest branch is much flatter than for the uniform gas. The explanation for this behaviour 
is once again found in Appendix B where we show that the long-wavelength frequency of 
plasmons propagating in the x direction is given by u 2 (q x ,q y = 0) = 2irn2D(qx/mx), where 
m x is the effective band mass at the zone center. Since m x increases with Vo, the dispersion 
of the plasmon becomes weaker with increasing modulation. 

In Fig. [| we present the g^-dispersion analogous to that shown in Fig. |5] for the case of 
weak modulation. As expected, the characteristic ^Jq^ behaviour is seen to persist (the "1-6" 
branch), but now the 2D plasmon crosses the lowest q x = ir/a mode (the "2-8" branch). 
At these crossing points, the dispersion of the plasmon band as a function of q x is weak 
since the frequencies at q x = and q x = 7t/a are coinci dentally equal. As for the case of 
weak modulation, we also see the lowest pair of q x = n/a modes (branches "2-8" and "3-7") 
crossing, implying that the lowest zone boundary gap closes at some finite value of q y . This 
again is possible due to the different parities of these modes ("2" and "3") in the middle 
panel of Fig. IUJ. The same behaviour is now seen for the next pair of q x = ir/a modes which 
also have opposite parity. The parity of the modes is of course preserved as a function of q y , 
but this is not the function of q x . The mode density "1" at q y = in Fig. |1(] evolves 

continuously into the mode density "2" as q x is increased from to n/a, and the parity of 
the mode remains unchanged. On the other hand, at q y = 3n/a, the mode density "6" 
evolves continuously into "7" , and thus a change in parity is observed as the mode disperses 
with respect to q x . This change in behaviour is associated with the crossing of the pair of 
q x = n/a branches "2-8" and "3-7". 

We again see in Fig. |9] the merging of zone boundary and zone center mode frequencies 
at large q y values. This behaviour is even clearer than in the case of weak modulation and 
the explanation is the same. The mode densities of the lowest two modes at q y = 37r/a are 
labeled "6" and "7" in Fig. ITDJ the density profiles of these two modes have a very similar 



shape near the cell boundaries where the equilibrium density is smallest, and only differ in 
the relative sign of the density fluctuation from one cell to the next. The weak interaction 
between the density fluctuations localized at adjacent cell boundaries accounts for the near 
degeneracy of these modes. The same can be said of the pair of density fluctuations labeled 
"8" and "9". 

The nature of the mode densities at the lowest zone center gap is shown in the lower-left 
panel of Fig. [TUJ. The mode labeled "4" is an odd-parity mode and contributes a peak to 
the power absorption near uo ~ 0.25 a.u. as shown by the curve labeled "(b)" in Fig. [7]. It 
is this mode which eventually evolves into the center of mass mode in the limit of isolated 
ID wires. There is an additional weak peak at u ~ 0.5 a.u. corresponding to the next 
odd-parity mode at the zone center. It is clear that increasing the modulation has led to a 
decrease in the amplitude of the Drude peak at u = and a transfer of oscillator strength 
to the other zone center modes. 



D. Strong Modulation 

The system is in the strong modulation regime when the phase point lies above the solid 
curve in Fig. 0. In particular, the phase point "(c)" corresponds to an equilibrium density 
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which is made up of an array of isolated quantum wires. Since there is no appreciable density 
overlap between adjacent unit cells, we expect a relatively flat g^-dispersion for the lower 
plasmon bands. In Fig. [11], we show the g x -dispersion for this case of strong modulation. 
Notice that the lowest plasmon branch has been pushed down to ui ~ (because of the large 
effective band mass m x ) and is dispersionless. Somewhat surprisingly, the next higher branch 
exhibits a dispersion which is larger than the next three branches lying between uj = 0.5 
and uj = 0.9 a.u. This would not be the expected behaviour in a tight-binding model of 
electronic band structure in which the dispersion of the bands increases with increasing 
energy (see Fig. ^). This conventional band structure effect is beginning to be evident at 
the top of Fig. |TT| As we shall see, the anomalous dispersion of the first excited band is due 
to inter-wire Coulomb interactions. 



The dispersion of the plasmons in the g^-direction is shown in Fig. [12] In contrast to 
the weak and moderately modulated cases, there are now two gapless modes present in the 
spectrum^. The long-wavelength ^Jq^ character of the q x = 2D plasmon "1" is still clearly 
apparent but there is a second mode "2" with q x = ir/a which has a linear dispersion!!. This 
branch originates from the "2-8" branch in Fig. ^ and in fact corresponds to a ID plasmon 
propagating along each wire, but with a density fluctuation which changes sign from one 
wire to the next. The two lowest curves in Fig. [12] define the plasmon band for an array 
of quantum wires and is qualitatively similar to the results found previously for a quantum 
wire superlattice based on the RPA&0 or hydrodynamic models!. 

The next mode "4" in Fig. [12] is the center of mass (CM) mode for a ID wire and the 
"4-9" and "3-7" branches define the CM plasmon band. This band is qualitatively similar 
to that found by Eliasson et al, Ref. [ |[. The density fluctuation at q x = and q y = is 
shown by "4" in Fig. [13] and, to an excellent approximation, is simply the derivative of the 
equilibrium density profile shown by curve "(c)" in Fig. II This is the expected behaviour 



of the 'sloshing' CM mode for parabolic confinement As one disperses along q x to "3" 
at q x = ir/a, or along q y to "9" at q y = 3ir/a, one sees that the density fluctuations at all 
positions are essentially the same. Thus, this mode is an intra-wire mode with very little 
interaction between the wires. 



Another interesting feature in Fig. [12| is the way in which the plasmon band merges with 
the CM band for large q y . This is reminiscent of the mode dispersions found in 3D parabolic 
wellst^l where the surface plasmon merges with the CM mode at higher wavevectors. The 
reason for the merging is the same as in the weak and moderate modulation cases; the density 
fluctuations shown by "7" and "8" in Fig. |T3] are symmetric and antisymmetric versions of 
each other with very similar density profiles at the edge of the wires. At this large q y value, 
the density fluctuations on either side of the wire interact with each other only weakly, as 
do the fluctuations on different wires. 



We turn next to the explanation of the finite bandwidth of the "3-4" branch in Fig. [Tl 
The origin of the dispersion is the dipole-dipole interaction between the wires. We consider 
the limiting case where the wire width W is small compared to the inter-wire separation a. 
For q y = 0, the electric field experienced by the n-th wire in the dipole approximation is 

E n = ^£'^ , (60) 
where p n is the dipole moment per unit length and the sum over s excludes the s = term. 
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We now suppose that the dipole moment has a plane wave modulation in the x-direction: 
Pn — Poe lQxna . In this case, the electric field is given by 

E n = • (61) 

a 2 g s 2 

The quantities p n and E n are connected through the relation p n = a{uj)E n where a{oj) is 
the dipole polarizability of the wire. Since the confining potential for the wire is roughly 
parabolic, the polarizability is given by a(u>) = A/(co> 2 — uj 2 ), where A = an-iD is the line 



density of the wire and uj q ~ y27r 2 Vo/a 2 is the frequency of the harmonic potential. Using 
this result in Eq. (|BT1), we find the following dispersion relation for the CM band: 



2 o 4A ^ cos(q x as) 



uj = uj n — 



l 2 i s 2 



(62) 



This result is entirely consistent with the observed q x dispersion for the "3-4" branch shown 
in Fig. 0. 

Finally, we discuss the FIR power absorption in Fig. [7[ Each successive curve corresponds 
to a modulation increase of AVo = 1.0 a.u., starting at Vq = 2.0 a.u. The curves labeled 
"(a)" and "(b)" were discussed earlier in the context of weak and moderate modulation, 
respectively. We note that an increase in the strength of the modulation transfers oscillator 
strength from the Drude peak to the higher q x = modes having a nonzero dipole moment. 
The two peaks which grow in intensity between curves "(a)" and "(b)" evolve from the 
q = 2ir/a and q = Arc /a 2D bulk plasmons. Both peaks are red-shifted with increasing 
modulation with respect to their positions in the uniform gas limit. Near the curve labeled 
"(b)", this trend is reversed and the peak frequencies begin to increase with further increase 
in modulation. At the same time, the oscillator strength of the lower peak continues to 
increase, while that of the higher frequency peak begins to decrease. Beyond this point we 
also see a rapid reduction in the strength of the Drude peak at uj = 0. 

The change in behaviour of the peak positions and intensities near the curve "(b)" occurs 
as the 2D- ID phase boundary in Fig. § is being crossed. Thus, a signature of the 2D — > 
ID transition is the observation of a minimum in the fundamental resonance frequency as a 
function of modulation, and the disappearance of the Drude peak. Beyond this minimum, 
the fundamental resonance grows in intensity as it evolves into the CM (Kohn) mode^TEil 
for parabolic confinement. For large modulations, the quantum wires are well-separated 



and sit in an effectively harmonic potential with frequency luq ~ y27r 2 Vo/a 2 , which explains 
the increasing trend observed in Fig. [7|. This also explains why the fundamental resonance 
eventually exhausts the total FIR oscillator strength. 

As an overview of our results, we plot in Fig. O the zone boundary and zone center 



frequencies at q y = as a function of the strength of the modulating potential. Only the 
lowest ZB modes exhibit a linear separation of the frequency with Vo and is in accord with 
our perturbation theory results. The lower branch is seen to tend to zero as Vq increases 
as this mode evolves into the ZB (or out-of-phase) plasmon of a ID wire superlattice. The 
lowest ZC mode (indicated by solid circles) tracks along the CM mode in Fig. ^ and, as 
discussed above, goes through a minimum at the transition from 2D to ID behaviour. This 
behaviour is qualitatively similar to that observed0@, although it should be noted that in 
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the experiments the average density of the gas decreased as a function of the gating potential, 
whereas our results are for a constant density. The CM mode crosses the upper branch of 
the second ZB mode at Vo/Ef ~ 20 and the pair then defines the CM band with its dipole 
induced dispersion. For the higher modes we see the merging of successive pairs of ZC and 
ZB modes in the ID limit. The increasing frequency of these modes is once again due to the 
increasing curvature of the confining potential. Fig. [14] extends the variation of the mode 
frequencies with Vo into the localized ID regime which is beyond that given previouslyB0. 
In particular, we do not find the unphysical result claimed in Ref. [ [10] that the system does 
not support plasmons in the strong localization limit. 



V. CONCLUSIONS 



In this paper, we have demonstrated that the Thomas-Fermi-Dirac-von Weizsacker ap- 
proximation provides a realistic description of the collective excitations in a modulated 
2DEG. We have presented a detailed investigation of the equilibrium properties, and have 
numerically mapped out the parameter space that defines the transition from 2D to ID 
behavior. We have also calculated the plasmon dispersions for propagation along and across 
the modulation direction in the weak, moderate and strong modulation limits. In agreement 
with earlier work on this problem, the modulation of the equilibrium density leads to the 
appearance of gaps in the plasmon dispersion in the direction of the modulation; as the 
modulation becomes stronger the bands become narrower and the gaps larger. For propa- 
gation in a direction perpendicular to the direction of modulation, we have shown that the 
long-wavelength 2D plasmon is unaffected by the modulation potential. However, at shorter 
wavelengths, the plasmon dispersions exhibit interesting behaviour as a function of q y and 
the explicit calculation of the mode densities provides a more complete understanding of the 
physical nature of these excitations. 

The power absorption in the long-wavelength limit has also been calculated for a range 
of modulation potentials. We have found that the oscillator strength is predominantly in 
the Drude peak for weak modulation and shifts to higher dipole modes as the modulation 
is increased. The reduction of the Drude peak is one useful indicator to gauge the effective 
dimensionality of the system. We have also shown that the 2D bulk plasmon at q — 2it/a 
evolves continuously with increasing modulation into the CM mode of a quantum wire 
superlattice. The frequency of the mode at first decreases, and then passes through a 
minimum as the 2DEG makes the transition from a continuous charge distribution to an 
array of isolated wires. This minimum is a second signature of the 2D — > ID crossover. At 
large modulations, the confining potential is effectively harmonic, and the FIR oscillator 
strength resides in the CM mode in accord with the generalized Kohn theorem. 

We are at present completing a detailed analysis of the magnetoplasma excitations in 
laterally modulated 2D systems. These results will be reported elsewhere. 
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APPENDIX A: CALCULATION OF PLASMON BAND GAPS 



In this appendix, we determine the eigenvalue spectrum of a laterally modulated 2DEG 
in the weak modulation regime. The eigenvalue problem that we wish to solve is 
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with 



! cg = B GG 'M G i G iic G " 

G'G" 

= N GG'C G i , 
G' 



N GG > = ^ B GG "M G n G i 

G" 



(Al) 



(A2) 



In the absence of an external modulation, the iV-matrix is diagonal. Denoting the matrices 
in this limit by a superscript 0, we have 



B GGI = ~{{q x + Gf + q 2 y )5 GG> 



Ml 



GG' 



Arm 



2D 



(q x + Gf + q\ 



+ ±C x n 2D - \C Z ^> + ^ ((& + G) 2 + q; : 



'GG' 



(A3) 



With these results, we have at q y — 0, 

n gg > = B GG „M G „ G , 



G" 



2nn 2 D\qx + G\ + [2Cin 2 D - -C 3 V^2~d) (Qx + G)' 
= ul{q x + G)5 GG i . 



>GG' 



(A4) 



This equation gives the plasmon frequency in Sec. |IV A| . With G taking on all possible 
values, q x can be restricted to the first Brillouin zone, —ixja < q x < ix/a. 

The introduction of an external modulation lifts the degeneracies in the uniform gas 
spectrum at the zone center (q x = 0) and at the zone boundary (q x = n/a). At the zone 
center (ZC), the reciprocal lattice vectors coupling degenerate modes are given by 



0, G\ = n—, G 2 
a 



7T 



-n- 



where n is an even integer, while at the zone boundary (ZB) 



q x = ~, G 1 = (n-1)-, G 2 
a a 



(n + 1) 



71 



(A5) 



(A6) 



where n is an odd integer. At either ZC or ZB, \q x + G\\ = \q x + G 2 \ = nir/a = q n and 
{G x - G 2 ) = Irmja = 2q n . 

We now apply degenerate perturbation theory to obtain first-order corrections to the 
frequencies at these points in the first Brillouin zone. Retaining only the two degenerate 
modes of interest, we have the pair of equations 
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u 2 c Gl = N GiGi c Gl + N GiG2 c G2 
u 2 c G2 = N G2Gi c Gl + N G2G2 c G2 



(A7) 



where N GiG2 = N GiGi is the lowest order off-diagonal correction to the iV-matrix. It is given 
by 



N hiG 2 = BgxGiMq^ + B GiG2 M G2G2 , 



(A8) 



where 



M GlG2 



q 



B, 



G1G2 



n —i/j [2q n ] + ^n [2q n ] - ^^ [2g„] ^=qlip [2q n ] 



Vn2D 



(A9) 



In the above, ip [2q n ] and no[2g n ] denote the n-th Fourier components of the groundstate 
wavefunction and charge density respectively Replacing ip [2q n ] by n [2q n ]/2y/n 2 D, Eq. (|A8[) 
can be cast in the form 



3 3 

-ul(q n ) + -C 3 V%d^ - ~\ w q* 



n [2q n ] 
n [2q ] 



(A10) 



This implies that the eigenvalues of (|A7| ) are given by 
ool(q n )=N°[q n }±\N 1 [q n }\ 

3 3 

= ^o(9n) ± ul(q n ) - -C 3 ^/n^>q* + -\ w q* 



n [2q n 



M 2 9o] 



(All) 



where A^°[g n ] = N G G and A^ 1 [g n ] = N GiG2 . The size of the gap at the ZC or ZB is given 
by 2|A? 1 [g n ]|. We have tested the validity of ( |A11| ), and in the regime where the pertur- 
bative results are valid, we find excellent agreement between the analytical and numerical 
calculations as illustrated in Fig. III. In the Hartree approximation, we obtain the simple 
relation 



w±(g„) = ul{q n ) 1 ± 



n [2q n ] 



™o[ 2 <?o] 



(A12) 



This result is entirely consistent with the one obtained by Krasheninnikov and ChaplikB in 
the long- wavelength limit. 



APPENDIX B: DETERMINATION OF LONG WAVELENGTH DISPERSIONS 

A persistent feature of the dispersion relations we have calculated (see Figs. f|, || and |TTD , 
is that the long-wavelength limit of the lowest plasmon has the characteristic y/q dependence 
of the 2D-bulk plasmon, regardless of the modulation. In this section, we show analytically 
that the long-wavelength behaviour of the lowest plasmon branch is ^fq for any modulation 
amplitude. 
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Once again, the eigenvalue problem that we are solving has the form 



j Cq — ^2 Bcg'Mg'G"Cg" 



G'G" 
1 



^ h-GG' Mg> G" C G" 



(Bl) 



where haG' is the Hamiltonian which, for q y = 0, determines the groundstate von Weizsacker 
wave function t/> I n the long-wavelength limit, the dominant contribution to the M- 
matrix comes from the G" = component of the Hartree matrix in Eq. (1311), 



lim M GG ' ^ lim Mq G , 



An 



(B2) 



q^o q^o — qa 

With this result, and noting that J2c ^GG'V'ot^"] = £ o( < l)i ) o[G], Eq. ([Bl]) can be written as 



47J- 

X w uj 2 c g = — e (ci)ij [G} YsMG")^" ■ 
qa q„ 



(B3) 



The eigenvalues u 2 of ( |B~3"| ) are thus given by 



2 Aim 2D i v 

u {qx,q y ) = — - — £o(q) 



(B4) 



since n 2 £i = a 1 X) G l^ot^] P- For small q we have s (q) ~ ^ + q%\, where m x is the 

effective band mass due to the confining potential in the modulation direction. Of interest 
here is the determination of the dispersions in both the q x = 0, q y — > and q y = 0, q x — > 
limits. In these two cases we have 



w 2 (<?x = 0, ~ 2'Kn 2 Dq y 



(B5) 



^ (ftc, ?y = 0) ~ — — q x 



m 7 



(B6) 



Eq. (|B5|) illustrates the fact that for any modulation, the long-wavelength behaviour of 
the 2D plasmon is ^fq~y~. Eq. ( (B~6|) on the other hand indicates that the ^fq^ dispersion is 
suppressed with increasing modulation as a result of the effective mass m x . 
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FIGURES 



FIG. 1. Equilibrium density profiles of Eq. (||) with w ex t given by Eq. @. The curves labeled 
by (a), (b), (c), and (d) correspond to Vq/Ej = 6.4, 19.1, 28.6 and 41.4, respectively. The mean 
density in all cases is t%2d = 0.1. 

FIG. 2. "Phase" -diagram for the 2D — ► ID transition. The solid and dashed curves include 
all interactions with 7%2D = 0.1 and U2D = 1-0 respectively. The open circles correspond to the 
noninteracting phase curve (i.e. a confined noninteracting Fermi gas) which is independent of the 
average 2D density. 

FIG. 3. The calculated TFDW energy bands along the direction of modulation. This figure is 
evaluated for a moderate modulation corresponding to "(b)" in Fig. |l[ 

FIG. 4. The dispersion of mode frequencies with wave vector q x for the weakly modulated 
2DEG. The solid curve is for the uniform 2DEG. A full explanation of the encircled numbers is 
given in the text. 

FIG. 5. The dispersion of mode frequencies with wave vector q y for the weakly modulated 
2DEG. Here q x = (solid circles), and q x = ir/a (open circles). A full explanation of the encircled 
numbers is given in the text. 

FIG. 6. The mode densities in the weak modulation regime for various bands. The left panels 
are evaluated at fixed q y = and the right panels at evaluated at q y = 3n/a. The sequence for the 
panels from top to bottom is q x = 0, q x = it /a, and q x = 0. See the text for full details. 

FIG. 7. Power absorption as a function of frequency. The curves are ordered from weakest 
(bottom) to strongest (top) modulation and are offset by AVb = 1-0 a.u. The labels (a), (b), (c) 
and (d) are consistent with the notation adopted in Figs. |l| and ||. 

FIG. 8. As in Fig. ||], but for moderate modulation. 

FIG. 9. As in Fig. but for moderate modulation. 

FIG. 10. As in Fig. |6|, but for moderate modulation. 

FIG. 11. As in Fig. |], but for strong modulation. 

FIG. 12. As in Fig. §, but for strong modulation. 
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FIG. 13. As in Fig. ||, but for strong modulation. 



FIG. 14. Zone boundary (open circles) and zone center (filled circles) frequencies (in a.u.) at 
q y = as a function of the strength of the modulating potential. The solid curves are the solutions 
given by Eq. ( |AlT|) for n = 1. 



FIG. 15. The first three Fourier coefficients of the equilibrium charge density as a function 
of modulation strength. The curves are normalized with respect to the average areal density 
Hid- Here G n = 2mr/a, with n = 1 corresponding to the figure inset (solid line), n = 2 the 
short-dash-long-dash curve and n = 3 the dashed curve. 
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